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Abstract 

We used a torsion pendulum containing rs 10 23 polarized electrons to search new interactions that 
couple to electron spin. We limit CP-violating interactions between the pendulum's electrons and 
unpolarized matter in the earth or the sun, test for rotation and boost-dependent preferred-frame 
effects using the earth's rotation and velocity with respect to the entire cosmos, and search for exotic 
velocity-dependent potentials between polarized electrons and unpolarized matter in the sun and 
moon. We find CP-violating parameters \gpgs\/{hc) < 9.4 x 10~ 37 and \g%g^\/(hc) < 1.2 x 10~ 56 
for A > 1AU. We test for preferred-frame interactions of the form V = —cr e A, V = —B<r e -v/c, 
or V = —J2 a iCijVj/c, where v is the velocity of the Earth with respect to the CMB restframe 
and i,j represent the equatorial inertial coordinates X,Y and Z. We constrain all 3 components 
of A, obtaining la upper limits |Ax,y| < 1.5 x 10~ 22 eV and \A Z \ < 4.4 x 1CT 21 eV that may 
be compared to the benchmark value m 2 /Mpi anc k = 2 x 1CT 17 eV. Interpreting our constraint 
on A in terms of non-commutative geometry, we obtain an upper bound of (355 /gut) 2 on the 
minimum observable area, where /gut = ?ic/(10 16 GeV) is the grand unification length. We find 
that \B\ < 1.2 x 1CT 19 eV. All 9 components of C are constrained at the 1CT 17 to 1CT 18 eV level. 
We determine 9 linear combinations of parameters of the Standard Model Extension; rotational- 
noninvariant and boost-noninvariant terms are limited at roughly the 10 -31 GeV and 10 -27 GeV 
levels, respectively. Finally, we find that the gravitational mass of an electron spinning toward the 
galactic center differs by less than about 1 part in 10 21 from an electron spinning in the opposite 
direction. As a byproduct of this work, the density of polarized electrons in Sm C05 was measured 
to be (4.19 ± 0.19) x 10 22 crn^ 3 at a field of 9.6 kG. 

PACS numbers: 11.30.Cp,12.20.Fv 
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I. INTRODUCTION 



This paper describes constraints on possible new spin-coupled interactions using a torsion 
pendulum containing «lx 10 23 polarized electrons. Several rather different considerations 
motivated this work. General relativity (the classical theory that forms the standard model 
of gravity) does not in itself make predictions about the gravitational properties of intrinsic 
spin-a quantum mechanical effect with no classical analog. Because the classic equivalence- 
principle and inverse-square-law experiments all used unpolarized test bodies they can shed 
no light on this issue. In addition, these "fifth force" experiments were completely insensitive 
to the purely spin-dependent forces arising from the first-order exchange of unnatural parity 
(0~, 1 + , etc) bosons. Lastly, the spin pendulum provided a means to search for a new class 
of preferred-frame effects that involve intrinsic spin. 

We established very tight bounds on long-range CP- violating interactions and on velocity- 
dependent forces by searching for the effects of unnatural-parity bosons exchanged between 
the electrons in our pendulum and unpolarized matter in the earth and in the sun. Similarly, 
we made sensitive searches for preferred-frame effects defined by the entire cosmos. In 
this case, we checked whether the spins in our pendulum preferred to orient themselves 
in a direction fixed in inertial space, or if they had a generalized helicity defined by their 
velocity with respect to the rest-frame of the cosmic microwave background. Finally, we 
investigated non-commutative space-time geometries and obtained constraints that lie far 
beyond the reach of any proposed accelerator. In each case, our bounds are interesting 
because of their high sensitivity. A Letter on the first results from this apparatus has 
already appearedfl]. This paper, which supercedes Ref. [1], reports additional data and more 
powerful constraints, provides a complete account of the experimental work, and presents a 
more extensive discussion of our constraints. 

A. Cosmic Preferred Frames 

1. General Considerations 

It seems natural to us, when testing possible Lorentz- violating preferred-frame scenarios, 
to let the cosmic microwave background (CMB) define the fundamental inertial frame. We 
consider, in turn, three types of preferred-frame effects involving spin. First, we consider an 
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interaction that violates rotational invariance and produces a laboratory potential 

V = -a-A, (1) 

where A is a preferred direction with respect to the cosmos as a whole Next, we consider 
a helicity-generating interaction that violates boost invariance as well, i.e. a lab-frame 
interaction 

V = -Bo- ■ v/c , (2) 

where v is the velocity of the spin with respect to the CMB rest-frame. This velocity has 
two significant components-the velocity of the earth with respect to the sun (v@/c ~ 10~ 4 ) 
and the velocity of the sun with respect to the CMB rest frame (i>cmb/c ~ 10~ 3 ); the small 
velocity due to the earth's rotation is neglected. 

Finally, following the general ideas of Kostelecky and coworkers[2], we consider a tensor 
helicity-generating term that produces a laboratory interaction 

^ = -£>i-<V (3) 

We take Cartesian equatorial coordinates for our inertial frame, where Z points North along 
the earth's spin axis, X points from the earth to the sun at the 2000 vernal equinox, and 
Y = Z xX. 

These possible spin-dependent cosmic preferred-frame effects differ from the known CMB 
preferred-frame effect in two important ways. First, if A is a time-even polar vector, B a 
time-even scalar and C a time-even tensor (untestable assumptions even in principle) , then 
the interactions in Eqs. 1, and either 2 or 3 violate time-reversal and parity, respectively. 
Second, if such spin-dependent preferred-frame effects exist, they could be observed within 
an opaque box that shielded out any normal information about the outside world. Although 
we naively expect that any such preferred-frame effects to be suppressed by the Planck scale 
and therefore of order m^/M-p^^ k 2 x 10~ 17 eV, our results are sufficiently sensitive to 
probe such tiny effects. 

2. The Standard Model Extension effective theory 

Kostelecky and coworkers [2] have developed a preferred-frame scenario in which spin-1 
and spin-2 fields were spontaneously generated in the early universe and were subsequently 
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inflated to enormous extents. These fields clearly violate Lorentz symmetry but in a very 
gentle way. Rotations and boosts of an observer remain Lorentz invariant, but the same 
operations on a particle (but not the Universe) obviously do not respect the invariance 
principle. The Lorentz non-invariance invalidates the Pauli-Luders theorem, allowing Kost- 
elecky et al. to construct a field theory with CPT-violating effects (the Standard-Model 
Extension or SME) that has been widely used to quantify the sensitivity of various CPT 
and preferred-frame tests involving photons, leptons, mesons and baryons (Ref. [3] has an 
extensive set of references to this work). 

Our A coefficient is identical to the SME coefficient b e which contains contributions from 
CPT-violating as well as CP- violating terms [4]. The SME has similar terms to our C 
coefficients, but they are conventionally expressed in sun-fixed coordinates. To facilitate 
SME analyses and to compare our results with other work based on SME analyses, we also 
quote results in which the sun's velocity with respect to the CMB rest frame is ignored. 

3. The Ghost Condensate dynamical theory 

Arkani-Hamed and colleagues have developed a consistent dynamical mechanism for the 
spontaneous breaking of Lorentz symmetry that, in effect, proposes a pervasive, massless, 
aether-like fluid consisting of Nambu-Goldstone bosons associated with the broken time- 
diffeomorphism symmetry[5, 6]. The bosons form a "ghost condensate" because the as- 
sumption of translational invariance demands a negative kinetic term in the Lagrangian. 
They form a condensate because higher-order terms are needed to stabilize the negative ki- 
netic term. This fluid has remarkable properties. Unlike a classical aether, it can be excited. 
It modifies Newtonian gravity at long ranges and late times in the evolution of the universe. 
It behaves like a cosmological constant and is capable of accelerating the expansion of the 
universe, but maximally violates the Equivalence Principle. It defines a preferred frame that 
converges with the rest frame of the CMB as it reaches equilibrium, and as it evolves to its 
final state it mimics dark matter. Cosmological tests of gravity limit the diffeomorphism 
symmetry-breaking scale M to be less than 10 MeV. Interesting physics is possible within 
this constraint; if M were ~ 1CT 3 eV the condensate could drive the observed acceleration 
of the universe; if M were ~ 1 eV it would mimic the observed dark matter. Particularly 
interesting in the present context, a fermion moving with respect to the ghost condensate 
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experiences an interaction 

M 2 

V=—vv/c, (4) 

where F is a mass scale associated with the coupling of the fermion to the condensate, and 
v is the fermion's velocity with respect to the condensate, which we assume to be at rest in 
the frame where the CMB is essentially isotropic. Equation 4 has the same form as Eq. 2. 

B. Non-commutative spacetime geometries 

Preferred- frame effects also occur in non-commutative space-time geometries [7, 8] that 
can arise in D-brane models. In these geometries, the space-time coordinates x M do not 
commute, but instead satisfy 

[x^Xj,] = , (5) 

where Ml/ is a real, antisymmetric tensor that has dimensions of length squared and |0| 
represents the smallest "observable" patch of area. It is often assumed that 0o« = to avoid 
problems with causality. We make the usual additional assumption that is constant over 
the space-time region spanned by our experiment, so that it defines a preferred direction 
rf = e' ljk <S)j k . One consequence of Eq. 5 is that a spinning electron experiences a CP- violating 
interaction [8, 9] 

£ cff = -mlf (—J e""^vV , (6) 

where A is a cutoff conventionally taken to be ~ 1 TeV, and a is the fine structure constant. 
Equation 6 is equivalent to 

^ = ^(l^v) 2 3 - 33 >< 1036 S 161 ' (7) 
which has the same general form as Eq. 1. 

C. Forces from exotic boson exchange 

The exchange of pseudoscalar particles produces spin- dependent forces that vanish be- 
tween unpolarized bodies. As a consequence, there are few experimental constraints on such 
forces. Moody and Wilczek[10] discussed the forces produced by the exchange of low-mass, 
spin-0 particles and pointed out that particles containing CP-violating J w = + and J w = CT 
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admixtures would produce a macroscopic, CP-violating "monopole-dipole" interaction be- 
tween a polarized electron and an unpolarized atom with mass and charge numbers A and 
Z 
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(8) 



where = h/(Xc) is the mass of the hypothetical spin-0 particle, gp and gs are its pseu- 
doscalar and scalar couplings, and g$ = Z{g% + gs) + (A — Z)g$. For simplicity, we assume 
below that gg = g$ = g^ and g$ = so that g$ = Ag$ ; constraints for other choices of the 
scalar couplings can be readily obtained by scaling our limits. 

Recently Dobrescu and Mocioiu[ll] classified the kinds of potentials, constrained only by 
rotational and translational invariance, that might arise from exchange of low-mass bosons. 
Our work is sensitive to 3 of their potentials; in addition to a potential equivalent to Eq. 8, 
we probe the two velocity-dependent potentials 
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where v is the relative velocity in units of c. Both terms can be generated by one-boson 
exchange in Lorentz-invariant theories. The parity-conserving f± term can arise from scalar 
or vector boson exchange, while the parity-violating f v term can be induced by spin-1 bosons 
that have both vector and axial couplings to electrons, giving /„ = 2g e A gy . We constrain 
the parameters of Eq. 9 by studying the interaction between the polarized electrons in our 
pendulum and unpolarized nucleons in the sun and moon. 



II. APPARATUS 



A. Rotating torsion balance 

This measurement used a substantially upgraded version of the Eot-Wash rotating torsion 
balance that was used for a previous test [12] of the Equivalence Principle. This device, in 
turn, was an improved version of the instrument described in detail in Ref . [13] . The essential 
features of this balance are shown in Fig. 1. Briefly, a vacuum vessel containing a torsion 
pendulum and its associated optical readout system was rotated uniformly about a vertical 
axis. The pendulum hung in a region of very low magnetic fields and gravity gradients, and 
the entire instrument was temperature-controlled (stability ±5 mK). Magnetic fields and 
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FIG. 1: Scale cross-section of the Eot-Wash rotating torsion balance; 1: torsion fiber, 2: auto- 
collimator, 3: Q21 gravity-gradient compensator, 4: nested magnetic and electrostatic shields, 5: 
turntable, 6: Q31 gravity-gradient compensator, 7: vacuum vessel, 8: spin pendulum, 9: ther- 
mal shield, 10: 3-axis Helmholtz coil system, 11: thermally controlled feet, 12: mode spoiler, 13: 
co-rotating slip-ring assembly, 14: upper fiber attachment mechanism. Two additional layers of 
magnetic shielding (not shown for clarity) are immediately inside and outside of the vacuum vessel. 

gradients were reduced by a stationary set of Helmholtz coils and 4 layers of co-rotating mu- 
metal shielding. Gravity gradients were canceled as described in Ref. [13], with a precision 
that was limited by the fluctuating water content of the ground outside our laboratory. 

For this measurement, we made several substantial improvements to the instrument used 
in Ref. [12]. The turntable was upgraded by a "feet-back" system, described in more detail 
in Appendix B, that kept its rotation axis vertical to better than 10 nanoradians. The 
autocollimator system that monitored the pendulum twist was upgraded from that described 
in Ref. [13] by using the "two bounce" geometry discussed in Ref. [14]. We also improved 
the resolution of the temperature monitoring system and upgraded the co-rotating magnetic 
shielding. 

The overall performance of the rotating balance is shown in Fig. 2. Except for signals 
at integer multiples of the turntable frequency caused by reproduceable irregularities in the 
turntable drive system, the noise in the twist signal is close to the thermal value expected 
from internal losses in the tungsten suspension fiber. Systematic errors associated with the 
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FIG. 2: Spectral power density in the twist signal. The horizontal scale shows the frequency in 
units of the turntable rotation frequency. The smooth curve shows the expected thermal noise 
density for this pendulum which had Q ~ 2000. At our signal frequency (shown by the arrow) 
the noise is predominantly thermal. Reproduceable irregularities in the turntable rotation rate 
produce peaks at integer values of ///tt- Spurious peaks at odd integer values were suppressed by 
combining data taken with 2 opposite orientations of the pendulum in the rotating apparatus. 

turntable signals were eliminated by combining data taken with different orientations of the 
pendulum within the rotating apparatus. 

B. Spin pendulum 

1. Design principles 

The heart of our apparatus is a spin pendulum, shown in Fig. 3, that contains a substantial 
number of polarized electrons while having a negligible external magnetic moment and high 
gravitational symmetry. The spin pendulum is constructed from 4 octagonal "pucks" . One 
half of each puck is Alnico 5 (a conventional, relatively "soft" ferromagnet in which the 
magnetic field is created almost entirely by electron spins) and the other half is grade 22 
Sm C05 (a "hard" rare-earth magnet in which the orbital magnetic moment of the electrons 
in the Sm 3+ ion nearly cancels their spin moment [15-1 7]). 

The trapezoidal elements of the pucks were fabricated by electric-discharge machining 
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and kept in precise alignment by aluminum frames. Thin plates glued to the sides of the 
Alnico pieces compensated for their lower density, p = 7.37 g/cm 3 versus 8.3 g/cm 3 for 
S111C05. After each puck was assembled, we magnetized the Alnico to the same degree 
as the Sm Co 5 by monitoring the external B fields as appropriate current pulses were sent 
through coils temporarily wound around the pucks. By stacking 4 such pucks as shown 
in Fig. 3, we placed the effective center of the spin dipole in the middle of the pendulum, 
reduced systematic magnetic-flux leakage, averaged out the small density differences between 
Alnico and SmCo 5 , and canceled any composition dipole that would have made us sensitive 
to violation of the weak Equivalence Principle. 

The pucks were surrounded by a gold-coated mu-metal shield that supported 4 mirrors 
equally spaced around the azimuth. The 107 g pendulum was suspended at the end of a 75 
cm long tungsten fiber which in turn hung from a mode spoiler that damped the pendulum's 
"wobble" , "swing" and "bounce" oscillations but had essentially no effect on the torsional 
mode. During the course of this experiment 4 different fibers were employed with diameters 
of either 28 jum or 30 /im. The upper attachment of the suspension fiber could be rotated 
to place any of the pendulum's 4 mirrors in the beam of the autocollimator system that 
monitored the pendulum's twist. 

We measured B , the field inside the pucks, by arranging 11 of the trapezoidal SmCo5 
elements in straight line, and used an induction coil to determine the field inside the Sm Co 5 . 
We found B = 9.6 ± 0.2 kG which was consistent with the supplier's[18] specification of 
9.85 kG. 

The absolute calibration of the twist to torque conversion was done by observing the 
pendulum's response to an abrupt change in the turntable rotation rate as shown in Fig. 4. 
The pendulum's free torsional oscillation frequency, fo = 5.379 mHz for this 30 /im diameter 
fiber, together with its calculated rotational inertia, determined the fiber's torsional constant 
k = 0.185 dyne-cm/radian. 

2. Estimating the spin content 

Polarized neutron elastic scattering studies[15] on SmCos have shown that the room- 
temperature magnetic moment of the Sm 3+ ion is extremely small, —0.04 fiB compared to 
the —8.97 hb moment of the five cobalts (see Appendix A). It is therefore an excellent 
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5.45 cm c- 

FIG. 3: (color online) Scale drawing of the spin pendulum. The light green and darker blue volumes 
are Alnico and SmCos, respectively. Upper left: top view of a single "puck"; the spin moment 
points to the right. Lower right: the assembled pendulum with the magnetic shield shown cut away 
to reveal the 4 pucks inside. Two of the 4 mirrors (light gold) used to monitor the pendulum twist 
are prominent. Arrows with filled heads show the relative densities and directions of the electron 
spins, open-headed arrows show the directions of B. The 8 tabs on the shield held small screws 
that we used to tune out the pendulum's residual (/21 gravitational moment. 

approximation that the magnetization of Sm C05 is due entirely to the Co. If we make the 
additional approximation (which we relax below) that the Co and Alnico magnetizations 
arise entirely from spins, then the net number of polarized spins in our pendulum is simply 

N p = ^Vrj (10) 

where B Q is the magnetic field inside a puck, R is the ratio of Sm to Co spin moments in 
room-temperature SmCo 5 (a negative quantity because the Sm and Co spins are directed 
oppositely), r\ = 0.65 accounts for the puck's octagonal shapes and V = 4.12 cm 3 is the 
total volume of SmCo 5 . As discussed in Appendix A, we used the Sm ion wavefunctions 
deduced from elastic neutron scattering[15] on SmCos to predict that at room-temperature 
the spin contribution to the Sm magnetic moment in SmCos is /is = +3.56/i#. Using data 
from this and other rare-earth/cobalt alloys, we estimate that the 5 Co's in SmCos have 
/is = —7.25fiB and obtain R = —0.49. 

On the other hand, a Compton scattering experiment on SmCos with circularly polar- 
ized synchrotron radiation [17], which directly probes the spin distribution rather than the 
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FIG. 4: Dynamical calibration of the torque scale, showing the pendulum twist as a function of 
time. At t = 377.7 s the turntable rotation frequency was abruptly changed from 0.55556 mHz to 
0.55611 mHz. The resulting change in the pendulum's oscillation amplitude and phase calibrated 
the angular deflection scale and determined the pendulum's free oscillation frequency /q. The 
smooth curve is the best fit. 

magnetic field, gives R = —0.23. We assume that R is equally likely to lie anywhere in the 
interval between —0.49 and —0.23 with a mean value and equivalent Gaussian spread of 

R = -0.36 ± 0.075 . (11) 

which implies N p = (0.79 ± 0.17) x 10 23 . 

However, as discussed in Appendix A, experiments have shown that the magnetization 
of Alnico and of Co in Sm C05 is not entirely due to spins as we assumed above. Defining / 
as the fractional contribution of spin to the total magnetic moment, we obtain an improved 
estimate 

N p = [f Co {\R\ - 1) + /Alnico] (12) 

= (0.96 ± 0.17) x 10 23 (13) 

where /amcc. — 0.953 ± 0.005 and /c = 0.81 ± 0.03 are taken from Eqs. A2 and A4 below. 

We show in Sec. VI F below that the Earth's rotation acting on the pendulum's angular 
momentum provided an internal calibration of N p in very good agreement with the value in 
Eq. 13. 
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III. EXPERIMENTAL SIGNATURES 



A. General form of the signatures 

The potentials of interest (Eqs. 1, 2, 3, 7, 8 and 9) all imply that the pendulum's energy- 
has the form 



where cr p represents the orientation of the pendulum's spin, and (3 has a value in the labo- 
ratory frame that may depend on time. This potential applies a torque r = N p <x p x (3 on 
the pendulum that can be detected by measuring its induced twist. However, only vertical 
torques, which are resisted by the very soft torsion spring constant, k, produce a measur- 
able twist; torques in the horizontal plane are resisted by the vastly stronger "gravitational 
spring" with an effective constant K grav = Mgs ~ 2.4 x 10 6 k where M is the pendulum mass 
and s = 2.72 cm is the vertical distance from the pendulum's center of mass to the fiber 
attachment point. As a result, the pendulum's spin is confined to the local horizontal plane 
and rotates along with the turntable about a vertical axis. Although the pendulum does 
not exactly follow the turntable angle because external torques on the pendulum twist the 
suspension fiber by angles 6 ^ 1CT 6 radians, to an excellent approximation, the pendulum 
twist, 9, as a function of turntable angle, 0, has the form 



Here (3± is the projection of (3 on the horizontal plane, 0o its azimuthal angle, and /3n and 
Pe are the North and East components of f3±. We use the convention that all angles are 
measured counterclockwise as seen from above, and that = points North. The procedure 
for converting our twist measurements into values for p N and Pe is discussed in Section V A 
below. 

The scenarios we test in this paper have signatures where (3 is fixed in the laboratory 
(Eq. 8), and others where (3 is fixed in inertial space (Eqs. 1 and 7) or depends on the 
orientation and/or velocity of the earth's motion around the sun (Eqs. 2, 3 and 9). These 
time- varying signatures were computed using astronomical formulae given by Meeus.[19] We 
neglected the earth's radius and the surface velocity associated with its spin in comparison 
to the astronomical unit and the earth's orbital velocity, respectively. 



E = -N p <t p • (3 , 



(14) 



6{<t>) = (P±N p /k) sin(0 o - 0) 

= -(A p //t)( / 9 N sin0 + /5 E cos0) . 



(15) 
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B. Lorentz-symmetry violation 



We constrained the rotational symmetry violating parameters of Eqs. 1 and 7 by searching 
for a f3 fixed in inertial space. Even though we are sensitive only to the two horizontal 
components of (3, the earth's rotation allowed us to probe all three components in the 
inertial frame, 

Ps(t) = [A x T^{t) + A Y T^{t)+A z T^{t)\ 

= -(A x cosQt + A Y sinQt)sm^ + A z cos$ 

(3 E (t) = [A x Tx(t) + A Y Ty{t) + A z T z (t)] 

= -A x sin Qt + A Y cos Qt . (16) 

The Txy Z (t) coefficients project the equatorial X, Y, Z system onto local N,E horizontal 
coordinates; \I> = 47.658 deg is the latitude of our laboratory, f2 = 27r/0.9972696d is the 
sidereal frequency and t is the local apparent sidereal time. 

We measured the boost parameters in Eqs. 2 and 3 using the motion of the earth with 
respect to the entire cosmos (neglecting the much smaller velocity of our laboratory with 
respect to the center of the earth). The velocity of the earth with respect to the rest frame 
of the cosmic microwave background (CMB) is 

V e = V Q + V C MB , (17) 

where v Q is the velocity of the earth with respect to the sun and v CM b is the velocity of the 
sun in the CMB rest frame. In equatorial coordinates 

VcMB = c(-1.190X + 0.253r-0.165Z)l(T 3 (18) 
v Q = VqX + w|(cos eY + sin eZ) , (19) 

where e = 23.45 degrees is the inclination of the ecliptic. The ecliptic coordinate z is 
perpendicular to the plane of the earth's orbit, x = X, and y = z x x so that 

Vq = —Rq cos A + XqRq sin A 

Vq = —Rq sin A — XqRq cos A q (20) 

where A (t), the true longitude of the sun, and -R©(t), the geocentric distance to the sun, were 
computed using expressions given by Meeus[19]. Converting v Q from ecliptic to equatorial 
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coordinates, we have 

/M0 = *E— ^(O, (21) 

i C 

and 

M) = Y.C^-^Tr(t). (22) 

Figure 5 shows the signatures of boost-symmetry violating terms with Cxx = 10~ 18 eV and 
with Cxz = 10~ 17 eV. Figure 6 shows how the /?e signature varies over the course of a year. 
The boost-violating signatures, which are modulated at the sidereal, annual and annual 
times sidereal frequencies, can be distinguished from a violation of rotational-invariance 
which have only a sidereal modulation. Note that the signatures of A z as well as Czx, Czy 
and Czz do not have daily modulations and must be extracted from lab-fixed /?n's. 



C. Exotic boson-exchange interactions 

1. Static potentials 

A CP- violating monopole-dipole interaction (Eq. 8) between the spin pendulum and the 
earth would give time-independent /3n and /3e values that depend upon the range, A, of the 
interaction. 

ge gN 

Pn,e = -j— -Tn,e(A) , (23) 



where 



WA) = e N , E .ft [(-L + -L) e-A] ^>£r , (24) 

where r points from the pendulum to a source mass element with density p, e N E is a unit 
vector pointing N or E, 

< 25 > 

and u is the atomic mass unit. This integral is, to within a constant factor, identical to that 
needed to evaluate the Eot-Wash Equivalence Principle results[13, 20], so we adopted the 
source integrations of those works. As explained in Refs. [13, 20], the integral can be reliably 
computed in two regimes: 1 m < A < 10 km (where information from local topography and 
boreholes is sufficient, and A > 1000 km where smooth, layered models of the entire earth 
are adequate. It is difficult to evaluate the integral in the gap region, and we have not done 
so. The integrals are displayed in Fig. 7. 
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FIG. 5: (color online) Daily signals produced by the Lorentz-boost violating term Cy. The daily 
signals can be essentially constant (i.e. with a weak annual modulation), modulated at the sidereal 
frequency, or modulated at both the sidereal and annual frequencies . Solid blue line: Cxx = 
10~ 18 eV; dashed red line Cyy = 10~ 18 eV. The points show the actual times at which data were 
acquired. 

2. Interactions between the spin pendulum and the sun 

A CP- violating monopole-dipole interaction (Eq. 8) between the spins in our pendulum 
and the sun would give modulated fa and fa values 

fa(t) = -{g e P g^/hc)I Q {\) cos a Q {t) cos 1& {t)/R 2 Q {t) 
M) = -(g e P g^/hc)I Q (X)cosa Q (t)sm lG (t)/R 2 Q (t) , 

(26) 
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FIG. 6: Expected /?e signal corresponding to Cxx = 10~ 18 eV. Each point is placed at the time 
of one of the cuts in our complete data set. Almost 4000 additional data points extending out to 
2825 on the horizontal axis are omitted to keep the figure readable. The annual modulation arises 
from v & , the steady component comes from vqmb- 

where a Q and 7 are the altitude and azimuth (measured Eastward from North) of the sun. 
For a uniformly dense sun of mass M & and radius r Q , the source integral becomes [21] 



KM, 
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(27) 



where £ = r /A. 

We searched for the long-range velocity-dependent boson-exchange interactions of Eq. 9 
using the earth's velocity around the sun. We needed, in addition to the velocity given in 
Eq. 20 and I Q (R, A), the term v x r — — \ Q R Q /c e z , and the source integral 
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8nu 



£ cosh £ — sinh £ 
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Rq' 


exp 


. A . 



(28) 



Then 



Ps(t) = f±I Q X& (cos a z cos 7 2 )/ (he 2 R & ) 

- f v I Q (v x cos a x cos ^ x + v y cos a y cos "fy)/R Q 
(fait) = f±I G \ G (cosa z sin>y z )/(hc 2 R & ) 

- f v I Q (v x cos a x sin ^ x + v y cos a y sin ^ y )/R Q 



(29) 
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FIG. 7: Earth source integrals denned in Eq. 24. The vertical axis shows /(A) divided by A. To 
obtain constraints on g%g% or gp§p, rather than on gf,gp , simply multiply /(A) by [Z//j]q « 0.48. 

where a(t) x ( y>z ) and r y{t) x (y,z) are the laboratory altitude and azimuth of the ecliptic unit 
vectors x, y and z. When evaluating the A-dependence of the couplings gpg$, f v and f± of 
Eqs. 8 and 9, we approximated the factors R in Eqs. 27-28 by their average value R — 1 
AU. 



3. Interactions between the spin pendulum and the moon 

We improved our /„ and f± constraints for A << 1 AU by using the moon's motion 
about the earth. The moon's coordinates and velocity were evaluated using expressions from 
Chapter 47 of Meeus[19] to compute generalizations of Eqs. 29 which are more complicated 
because the moon's motion does not lie in the ecliptic plane and is strongly perturbed by 
the gradients in the sun's field. 
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D. Gyro-compass effect 



Because the pendulum's magnetic flux was confined entirely within the pucks, the sta- 
tionary pendulum necessarily contained a net total angular momentum J = —S, where 
S = N p h/2 is the pendulum's net spin. The magnetic fields in the two materials, 

B z (SmCo 5 )/ii B = g s S z (SmCo 5 ) + 9l L z (SmCo 5 ) (30) 
B z (A\mco)/nB = g s S z (Ahdco) + ^L 2 (Alnico) , 

were equal and opposite so that the net magnetic field in the pendulum's z direction vanished. 
This requires that 

2(ST) + (Lf ) = , (31) 
where we make the excellent approximation that g s = 2. This implies 

(Jf) = -<5J*) • (32) 

To keep J fixed in the frame of the rotating earth, the fiber has to apply a steady torque, 
T, given by 

T-h = |n ffi x J h\ , (33) 

where £1$ is sidereal rotation frequency and n is the local vertical. The pendulum twists 
until kQ = —T • ft, i.e. J tries to point toward true North with a torque 

fl^lff, COS QQ . 

T • h = N p — = N p x (2.590 x 10" 39 Nm) , (34) 

where ^ is defined in Eq. 16. This is equivalent to a small negative (because J = —S) value 

= -1.616 x 10- 20 eV , (35) 

which formed a steady, precisely known background in our study of lab-fixed signals, and 
allowed us to calibrate the spin content of our pendulum. 

IV. MEASUREMENTS 

The data for this paper were taken in two distinct sets. Runs in set I were taken over 
a span of 18 months between August 2004 and January 2006. Data were taken with the 
autocollimator beam reflected by each of the pendulum's 4 mirrors. During this period, along 

19 



with acquiring normal spin-pendulum data, we made extensive studies of systematic effects 
from gravity-gradient and magnetic couplings and from turntable imperfections. Despite 
considerable effort, that included changing suspension fibers and mode spoilers, we were 
unable to understand a slowly- varying systematic effect, and as a result our laboratory-fixed 
constraints were dominated by poorly understood systematics. This work formed the basis 
of our brief publication [1]. 

We then decided to tackle the problem of slowly varying systematics by constructing the 
93 g "zero-moment" pendulum shown in Fig. 8. This aluminum pendulum had no net spin, 
was nonmagnetic and had the same moment of inertia as the spin pendulum. It coupled 
negligibly to gravity gradients because it had cylindrical symmetry and very small q2o, (feo 
and mass moments. In addition, we modified the fiber attachment to the pendulum so 
that we could reproduceably rotate the pendulum with respect to its suspension system. 
As discussed below, this device revealed an unexpected systematic effect associated with 
the suspension system that could be canceled by combining results taken with opposite 
orientations of the pendulum with respect to the suspension system. 

We added the new rotatable fiber attachment mechanism to the spin pendulum, and set 
II runs occured between September 2006 and March 2008. Data were taken with the light 
beam reflecting from each of the 4 mirrors and for 4 equally spaced angles between the 
pendulum and the suspension system. The results from set II gave much tighter lab-fixed 
constraints, and when combined with the set I data improved the sidereal and solar signals 
published in Ref. [1]. 

V. DATA ANALYSIS 

A. Extracting /3n and /3e from the pendulum twist 

Over the course of this experiment the turntable rotation frequency, /, was set at values 
between 3/ /29 and 3/o/20, where f is the free-oscillation frequency of the pendulum. We 
recorded the pendulum twist angle 9 as a function of the angle of the turntable (measured 
counter clock- wise from North), and converted it to torque, as described in Ref. [13]. Briefly, 
data were divided into "cuts" each of which contained exactly 2 revolutions of the turntable, 
lasted no more than 3800 s and typically contained 254 data points. A simple digital filter 
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FIG. 8: Scale drawing of the "zero-moment" pendulum. The pendulum body was made from 7075 
aluminum. The 4 vertical screws (shaded) were adjusted to trim away the residual 521 moments of 
the pendulum. The "ball-cone" device on top of the pendulum allowed us to rotate the pendulum 
with respect to the torsion fiber. The entire pendulum, including the mirrors, was coated with 
gold. 

was applied to suppress the free oscillations. The twist data from each cut were then fitted 
by 2nd-order polynomial "drift" terms (typical drift was 0.5 /xrad/h) and harmonic terms 

8 

9 ((f)) = (a£j cos ncj) — sin ncj)) . (36) 

n=l 

Our signal was contained in the n = 1 coefficients that we denote below as <2N(0d,0s) and 
a E(0d, 4>s), where cp d is the angle of the pendulum (normally its spin dipole) with respect to 
the turntable and cj> s is the angle between the pendulum and its suspension system. 
We converted the and cze signals for each cut into /3n and (3e values using 
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The fit parameters (djss, Se) are the twists associated with the pendulum itself, so that 

Pn,e = K^ N)E /iVp , (38) 
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FIG. 9: Data from a set of runs at (fid = 112.5 deg that spanned a duration of 113 h. The sidereal 
average of the signals has been set to zero. The dashed curves show the signal from a hypothetical 
A = (2.5 x 10 -20 eV)x which gives out-of-phase sine waves in (3n and (3e- The solid curves show 
the best sidereal fit, which yields A x = (-0.20 ±0.76) x 10~ 21 eV and A y = (-0.23 ±0.76) x 10~ 21 
eV. (Note that the hypothetical signal in Fig. 2 of Ref. [1] was actually A = (2.5 x 10~ 20 eV)x.) 

while (<5jvj, <5|) and (6^, 6%) are spurious twists associated with the orientations of the sus- 
pension system and turntable rotation-rate irregularities, respectively. For data taken with 
the 30/im fiber, <5 N E = 1 nrad corresponds to /?n,e = 1-2 x 10~ 21 eV. 

B. Astronomically modulated signals 

When analyzing our results for astronomically modulated signals, it was sufficient to 
ignore the last matrix product term on the right hand side of Eq. 37. We adopted the 
strategy of Ref. [13] and assumed that (3's modulated at daily, sidereal or yearly periods 
could be treated as essentially constant during any one cut. This approximation attenuated 
sidereal signals by at most 0.3% and was neglected. We used linear regression to fit these 
measured /3's with the expected signatures at the midpoint of each cut. Our input consisted 
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of 93 data sets (with durations ranging from 1 to 5 days) containing 8471 cuts and spanned 
a total of 1110 days. We minimized our sensitivity to slow drifts in our signals by zeroing the 
average values of the twist signals in each data set as well as any steady average value of the 
predicted signatures. The individual data sets were weighted by the inverse-squares of the 
spread of the /3's in that set. We constrained A in Eq. 1 by fitting our data with the two free 
parameters A x and A y given in Eq. 16. Figure 9 shows a typical data set and the fit. The 
helicity term in Eq. 2 was constrained by fitting Eq. 21 with the single free parameter B. An 
analysis with 6 free parameters provided the constraints on the 6 modulated boost-symmetry 
violating C terms in Eqs. 3 and 22 (only 6 of the 9 C^-'s produce modulated signals because 
spins pointing along z give a vanishing E signal and a steady N signal). Limits on the boson 
exchange potentials defined in Eqs. 8 and 9 obtained using the Sun as a source were found by 
fitting Eqs. 26 and 29 with the 3 free parameters g^g^, f± and f v . Because of the inclination 
of the earth's rotation axis, these torques have components modulated with a 24 hour period 
as well as annual modulations. In every case, the resulting fit had xV^ — 0-99. Because the 
individual data sets were normalized to y?jv = 1, a x 2 / u — 0.99 for the combined 93 data 
sets showed that there was no evidence for any anomalous behavior. The sidereal or solar 
modulation of the signatures eliminated many systematic effects that were fixed in the lab; 
as a consequence our bounds on Ax and Ay and on the 6 modulated C terms are tighter 
than those on A z or Czx, Czy and Czz which are based on the lab-fixed limits discussed 
below (see Tables VI and VII in Sec. VII). 

C. Laboratory-fixed signals 

1. Data set I 

We analyzed our /3n and /3e data for torques fixed in the lab frame by combining the signals 
observed for 4 equally-spaced angles, <pd, of the spin dipole within the rotating apparatus. As 
described in Ref. [1], we observed a scatter that was larger than our statistical uncertainties. 
This prompted us to study lab-fixed systematic effects with the "zero-moment" pendulum. 
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FIG. 10: (color online) Systematic error study using the "zero- moment" pendulum. Data were 
taken for 4 orthogonal directions of the pendulum with respect to the turntable and for two opposite 
orientations of the pendulum with respect to its suspension system and fitted with Eq. 37. The 
ellipses show the values and la uncertainties in the extracted signals associated with the pendulum, 
the turntable and the suspension system, respectively. Systematic effects associated with the 
pendulum itself were negligible. 

2. Systematic investigations with the zero-moment pendulum 

We analyzed data taken on all 4 mirrors and with two opposite orientations of the sus- 
pension system with respect to the pendulum using Eq. 37, where the 0d was denned by 
a ficticious zero mark on the pendulum. The results, shown in Fig. 10, revealed an unex- 
pected systematic effect associated with the pendulum suspension system, along with the 
known turntable systematic. Figure 10 shows that 5^ and 5| changed slightly every time 
we stopped the turntable and opened the vacuum system to change <p s . Therefore those 
parameters are shown as functions of <j) s in Eq. 37. There was no resolvable systematic effect 
associated with the pendulum. We therefore base all our lab-fixed constraints on the set II 
data, which did not suffer from the suspension systematic. 
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FIG. 11: (color online) Extraction of lab-fixed signals from the October 2006 data, showing the 
N and E components of the twist signals, obtained using Eq. 37. Ellipses indicate la statistical 
errors. The turntable systematic changed each time we opened up the vacuum system to change 
0dsi while the suspension systematic remained constant at a value very close to that seen in the 
"zero-moment" data plotted in Fig. 10. The red vertical lines shows the expected gyrocompass 
signal; the length of the line reflects the uncertainty in the R value. 

3. Data set II 



Spin-pendulum data were acquired with 4 orientations of the pendulum with respect to 
the rotating apparatus, and for 2 or 4 orthogonal orientations of of the suspension sys- 
tem with respect to the pendulum. Three such measurements were taken: Measurement 1 
(centered around October 15, 2006), Measurement 2 (centered around June 15, 2007), and 
Measurement 3 (centered around March 8, 2008). We used Eqs. 37 and 38 to extract from 
each of these measurement "true" fa and fa values that depended only on the orientation 
of the spin dipole. The extracted values of fa and fa for one of these measurements are 
plotted in Fig. 11. A highly significant signal was observed (see Table V below) that points 
within 0.2 ± 1.6 degrees of due South; it agrees well with the expected gyro-compass effects 
discussed in Sec. HID. 



25 



VI. SYSTEMATIC ERRORS 



The interactions in Eqs. 1—9 all produce lev twist signals. Any other mechanism that can 
produce a Iuj signal is a potential source of systematic error. We present below systematic 
error analyses for two types of such effects: lab-fixed signals whose amplitude and phase 
are constant in time and astronomical signals whose amplitude and phase vary as the earth 
rotates relative to the astronomical source. (Average daily variations of some important 
environmental parameters are given in Table I). In addition to such "false effects", we 
evaluate scale factor uncertainties in N p , the number of polarized electrons in our pendulum. 



TABLE I: Average daily components the systematic effects 



effect daily variation 

ambient magnetic field 2.0 ± 0.4 mG 

tilt 4.7 ±2.1 mad 

Q22 gravity gradient (1.10 ± 0.27) x 10~ 3 g/cm 3 

Q21 gravity gradient (5.0 ± 3.4) x 10" 4 g/cm 3 

temperature drift 0.34 ± 0.21 mK 

Iuj temperature 2.2 ± 1.7 [xK 



A. Turntable Rotation Rate 

We observed an 80 nrad leu signal that was very reproducible from cycle to cycle and 
independent of the orientation of the spin pendulum within the rotating apparatus. This 
arose primarily from a non-constant rotation rate caused by imperfections in the turntable 
bearings and appeared as a twist signal because we measured the angular position of the 
"inertial" pendulum relative to the rotating vacuum vessel. Although a servo system locked 
the output of the turntable's rotary encoder to a crystal oscillator to produce a nominally 
constant rotation rate, the finite gain of the servo loop and imperfections in the rotary 
encoder allowed residual twist signals to appear at the harmonics of the turntable rotation 
frequency as illustrated in Fig. 2. 
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We model the true phase, 0, of the turntable as 

<j> = wt+Y,<j> n e iu » t 

n=l 



where oj is the nominally constant angular frequency of the turntable, u n = nu, and <p n is 
the amplitude of the phase variation due to turntable imperfection. It is easy to show that 
(j) n will produce an angular twist signal, 9 n , at frequency u n (ignoring damping) given by 

On = 0n^-^ (39) 

where u p is the free torsional angular frequency of the pendulum. 



1. Astronomical signals 

We searched for time variations, A0„, of the turntable imperfections whose amplitudes 
and phases project onto the time dependence of our astronomical sources by fitting the time 
sequence of the 6* 4 and 9$ values to the basis functions of each astronomical signal. We chose 
#4 and 6*8 because signals at 4u and 8u do not change as the pendulum is rotated to each 
of its four positions within the vacuum vessel. We found no evidence for a modulation of 
9^ or # 8 at any of the astronomical frequencies. Furthermore, had there been a coherent 
modulation, by taking data at all four orientations of the spin pendulum within the vacuum 
vessel, the lu signal from A<pi would have been reduced by a factor of 0.22 (more data was 
taken on mirror 1 than on the other mirrors). Assuming that Acf>i ~ A0 4 ps A0 8 a s would 
be the case for a stiff point on the turntable bearing, and using Eq. 39 with the factor of 
0.22, the systematic error from A0i is less than 0.5% of the statistical error for any of the 
astronomical sources. This error is significantly smaller than the other systematic errors in 
Tables II and III and was dropped from further consideration. 



2. Lab-fixed signals 

We distinguished lab-fixed sources that couple to the spin dipole of the pendulum from 
leu signals that are independent of the spin dipole by taking data at four equally spaced 
orientations of the spin dipole within the apparatus. A systematic error can arise if 4>i in 
Eq. 39 changes as the orientation of the spin dipole is rotated. We searched for such a 
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TABLE II: Systematic and random errors in the astronomically modulated A, B and C signals 
evaluated in the CMB rest frame. Corrections were only applied for magnetic effects; the remaining 
effects were assigned systematic errors equal to the absolute value of the effects plus their la 
uncertainties. The total systematic errors are quadratic sums of the uncertainty in the magnetic 
correction and the systematics of the remaining effects. A 2.3% scale-factor uncertainty is not yet 

included. 

signal systematic error [eV] random error [eV] 



magnetic tilt grav. gradient temp, drift lu temp, total 



A x x 10 22 


+0.20 ±0.07 ±0.31 


±0.07 


±0.04 


±0.05 


±0.33 


±1.31 


A Y x 10 22 


-0.10 ±0.03 ±0.24 


±0.20 


±0.03 


±0.02 


±0.32 


±1.32 


B x 10 20 


-1.40 ±0.46 ±2.5 


±0.38 


±0.23 


±0.37 


±2.61 


±11.0 


Cxx x 10 18 


+0.95 + 0.31 ±0.43 


±0.17 


±0.08 


±0.05 


±0.57 


±2.16 


C X Y x 10 18 +2.29 ± 0.76 ±0.84 


±0.25 


±0.21 


±0.06 


±1.18 


±4.24 


Cxz x 10 18 


-3.61 ±1.19 ±2.58 


±0.99 


±0.31 


±0.22 


±3.03 


±10.16 


Cyx x 10 18 


-0.73 ±0.24 ±0.44 


±0.33 


±0.04 


±0.06 


±0.60 


±2.17 


C YY x 10 18 


-1.04 ±0.34 ±0.76 


±0.67 


±0.15 


±0.07 


±1.08 


±4.25 


C YZ x 10 18 


+4.02 + 1.33 ±2.46 


±1.60 


±0.21 


±0.32 


±3.24 


±10.23 



dependence by comparing the # 4 and 9g signals for the four orientations of the spin dipole 
and found no variation beyond the statistical error. Assuming that A0 X A0 4 A0 8 , 
Eq. 39 predicts that A<f>i produces a lab-fixed systematic of 0.11 x 10" 22 eV as listed in 
Table IV. 

B. Turntable Tilt 

If the rotation axis of the turntable is tilted from local vertical, the upper fiber attachment 
point flexes as the apparatus rotates, leading to a "tilt-twist" lu angular deflection of the 
pendulum. Imperfections in the mirrors of the system that measures our twist signal produce 
a related effect as the tilt causes the light spot to move on the mirrors. The tilt of our 
laboratory floor changes by a few /irad per day which, if not corrected, would produce 
a spurious leu signal that varies by up to 40 mad per day. The "feetback" control loop 
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TABLE III: Systematic and random errors in the heliocentrically modulated SME coefficients 
and one-boson exchange interactions. Corrections were only applied for magnetic effects; the 
remaining effects were assigned systematic errors equal to the absolute value of effects plus their 
la uncertainty The total errors are quadratic sums of the uncertainty in the magnetic corrections 
and the systematics of the other effects. A 2.3% scale-factor uncertainty is not yet included. 

signal systematic error [eV] random error [eV] 

magnetic tilt grav. gradient temp, drift Iuj temp, total 



A x x I0 22 


-0.05 ±0.02 


±0.30 


±0.06 


±0.05 


±0.05 


±0.31 


±1.44 


A Y x 10 22 


+0.08 ±0.03 


±0.21 


±0.17 


±0.03 


±0.02 


±0.27 


±1.44 


C' xx x 10 18 


±0.94 ±0.31 


±0.44 


±0.18 


±0.08 


±0.05 


±0.58 


±2.16 


C' XY x 10 18 


±0.73 ±0.24 


±0.64 


±0.23 


±0.12 


±0.06 


±0.73 


±1.98 


C'yx x 10 18 


-0.73 ±0.24 


±0.44 


±0.33 


±0.04 


±0.06 


±0.60 


±2.17 


C'yy X 10 18 


±0.70 ±0.23 


±0.56 


±0.08 


±0.16 


±0.11 


±0.64 


±1.99 


f ± /(hc) x I0 32 


-0.15 ±0.05 


±0.04 


±0.31 


±0.11 


±0.04 


±0.34 


±2.09 


f v /(hc) x 10 56 


-1.29 ±0.43 


±0.68 


±0.14 


±0.07 


±0.12 


±0.83 


±2.32 


&gg/(hc) x 10 37 


-0.10 ±0.03 


±1.60 


±0.40 


±0.31 


±0.29 


±1.70 


±8.53 



described in Appendix B aligns the rotation axis to within 10 nrad of the local vertical. 

We measure the tilt-twist feedthrough (the ratio of the spurious signal to the tilt) by 
applying known tilts to the apparatus; we find that the feedthrough is typically 3% with an 
amplitude and phase that depend upon the angular orientation of the pendulum within the 
apparatus. A conservative upper limit to the component of the tilt-twist feedthrough that 
mimics a coupling to the spin dipole of the pendulum (i.e. that rotates along with the spin 
dipole) is 4%. 

1. Astronomical signals 

We extracted any tilt component that mimicked an astronomical signal by fitting the 
outputs from the rotating tilt sensors to the basis functions of each of the astronomical 
signals. The largest projection, just over 3 standard deviations, was found for the Cxy 
basis function. We multiplied these tilt projections by the tilt-twist feedthrough, 4%, to 
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TABLE IV: Error budget for lab-fixed signals. A scale-factor uncertainty equaling 18% of the 
central value must be folded into the quadratic sum of the random and statistical errors. 



date 


10/2006 


uncertainty in f3 x 10 22 [eV] 
6/2007 


3/2008 


systematic effect 








tilt 


0.52 


0.88 


0.36 


PTaVlt.V PTfl H lPTlt, 

IA> V i V y £-) ^ CILL1L11 Li 


0.77 


0.56 


0.65 


temp drift 


0.14 


0.14 


0.20 


Iuj temp 


0.14 


0.14 


0.20 


turntable speed 


0.11 


0.11 


0.11 


magnetic 


1.06 


1.06 


0.51 


total systematic error 


1.43 


1.50 


0.95 


random error 


7.5 


5.6 


7.7 



determine the systematic error due to tilt. The results, given in Table II, show that control 
of the tilt-twist feedthrough is the largest source of systematic error for the astronomical 
signals. 

2. Lab-fixed signals 

The average tilt of the apparatus was 2.2, 3.7 and 1.5 mad for lab-fixed Measurements 
1, 2 and 3, respectively. The component of the tilt-twist feedthrough that rotated with 
the orientation of the spin dipole was no more than 2% for these data sets. The lab-fixed 
systematic error from tilt is then 2% of measured tilt, given in Table IV. 

C. Magnetic Coupling 

The spin pendulum had a residual magnetic dipole moment of 7 x 1CT 3 erg/G and was 
located within four layers of magnetic shielding. When the currents to the Helmholtz coils 
outside of the apparatus were reversed, a horizontal field of 374 mG was applied to the 
apparatus which produced a lu signal with a typical magnitude of 20 ±3 mad (f3 = 2 x 1CT 20 
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eV) . The spin pendulum was insensitive to reversal of the current in the vertical Helmholtz 
coil with a feedthrough of less than 2 ± 3 nrad for a 470 mG change in vertical field. 

1. Astronomical signals 

Before taking the data presented in this paper, as well as during the last 20% of the data, 
a triple-axis fluxgate magnetometer was placed within the Helmholtz coils (but outside of 
the heat shield) to monitor the stability of the magnetic environment. We were surprised to 
find an average daily variation of the horizontal magnetic field of 2.0 ±0.4 mG (see Fig. 12). 
Because the astronomical signals have a large daily component, the magnetic coupling to 
daily field variations could lead to systematic errors as large as 56% of the statistical error. 
To reduce this error, we corrected our data for the coupling to the daily variation of the 
magnetic field. No other corrections to the astronomical signals were required. 
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FIG. 12: Daily variation of the ambient magnetic field measured with a flux-gate magnetometer. 

The amplitude and phase of the magnetic coupling to the spin pendulum were found by 
reversing the Helmholtz coil currents, while the fluxgate measurements provided us with the 
amplitude and phase of the daily variation of the horizontal magnetic field. The calculated 
magnetic coupling was subtracted from each data cut before fitting the data to the astro- 
nomical signal basis functions. We assigned a systematic error for this procedure as follows. 
For each astronomical signal, we fitted the data with and without the magnetic correction 
to determine the magnitude of the correction (always less than 56% of the statistical error). 
There was a 15% uncertainty in the magnetic feedthrough calibration and a 29% uncertainty 
in the amplitude of the daily magnetic field variation (obtained from the difference between 
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the measurements before and at the end of the data sequence). We added these uncertainties 
in quadrature to obtain magnetic systematic errors equal to 33% of the magnetic corrections. 
Results are listed in Tables II and III. 

2. Lab-fixed signals 

An applied horizontal field of 374 mG gave a magnetic leu signal, Sm, of 22.2 ±1.4 nrad in 
Measurements 1 and 2, and 10.8 ±2.4 nrad in Measurement 3. The Helmholtz coil currents 
were adjusted to cancel the local field to within 0.5 mG with an uncertainty of 1 mG due 
to the daily variation of the horizontal field. We compute the magnetic lab-fixed systematic 
error from 1.5/374 x S M , i-e. 0.089 nrad (f3 = 1.06 x 10~ 22 eV) in Measurements 1 and 2 
and 0.043 nrad in Measurement 3. 

D. Gravity Gradients 

We use a spherical multipole basis[13] to characterize the mass moments, qi m , of the spin 
pendulum. The coupling of qu mass moments to Q n gravity gradient fields produce luu 
angular deflections of the pendulum that have the same dependence upon the orientation 
of the pendulum as would a coupling to the spin dipole. Although the spin pendulum was 
designed to have vanishing qu moments, imperfections produced small residual moments. 
The only significant gravity gradient coupling for this work was q 2 \Q 2 \ ■ This coupling was 
minimized by using a special gradiometer pendulum with a large g 2 i moment to measure 
ambient Q21 field. We mounted machined Pb/brass masses on a turntable just outside of 
the vacuum vessel to cancel the local Q 2 \ field by a factor of 200. Finally, we temporarily 
rotated the Pb/brass masses to double the gradient, measured the residual g 2 i moment of 
the pendulum, and adjusted tuning screws on the pendulum (see Fig. 3) to minimize the 
residual g 2 i mass moment. By canceling the Q 2 i field and minimizing the q 2 ± moment of the 
pendulum, the q 2 iQ2i coupling was reduced to tolerable levels. 
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1. Astronomical signals 



We were concerned that daily variations of the local Q21 held might couple to the residual 
q 2 i moment of the pendulum to produce a daily modulated leu signal that would project onto 
the astronomical signal basis functions. The Q22 held from local sources is typically larger 
than the Q21 field (because Q22 is maximum at a polar angle of 90° where Q21 vanishes) and 
falls off with distance with the same radial dependence as the Q21 field. The pendulum had a 
residual ^22 moment (uncompensated) that was 7 times larger than the residual 521 moment. 
We used the ^22^22 coupling to monitor the time variation of the Q22 field which for local 
sources will be correlated with changes in the Q21 field. The ^22^22 coupling produces a 2u 
signal. We fitted the time sequence of the measured 2uj signals to the astronomical basis 
functions to determine the projection of any Q22 variations onto those functions. We then 
ran a correlation analysis between the leu and 2uj signals over the entire data set to extract 
the linear slope between the leu and 2oj signals. That slope was found to be less than 3%. 
The gravity gradient systematic error in Table II is the projection of the 2eu signal onto the 
basis functions times 0.03. 

An independent analysis of the daily variation of the Q 2 i field was made using 8 day-long 
data runs with the gradiometer pendulum. The gradiometer pendulum had a q^i moment 
90 times larger than the residual moment of the spin pendulum. No evidence for a daily 
variation in the Q21 field was found with an uncertainty 4 times larger that obtained from 
the 522^22 analysis. 

2. Lab-fixed signals 

When the Pb/brass Q21 compensation masses were rotated to add to the local Q21 gra- 
dient, a change in loo signal of 70 nrad typically was measured. This coupling would change 
each time the ball-cone attachment on the spin pendulum was rotated because of small 
changes in the tilt of the pendulum relative to the ball-cone. The residual leu signal when 
the Pb/brass compensators canceled the local gradient was therefore less than 0.2 nrad, 
typically 20% of the statistical error for each ball-cone configuration. Because we knew the 
amplitude and phase of the total Q21 field and the amplitude and phase of the residual #21 
moment for each ball-cone configuration, we corrected the raw data for the residual q2iQ2i 
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coupling. This was the only correction applied to the lab-fixed signals and never exceeded 
30% of the statistical error for any configuration. The largest uncertainty in this procedure 
was the assumption that the total Q 2 i field at the pendulum did not change appreciably 
with time. We assign a systematic error to this procedure by taking the difference between 
the lab-fixed signals with and without the (721Q21 correction applied, and multiplying this 
difference by 50% to account for possible changes in the Q21 field. 

E. Thermal Effects 

The equilibrium twist angle of the spin pendulum is a strong function of temperature. 
Because these and other thermal effects are difficult to model, we took special care to 
stabilize and monitor the thermal environment of the apparatus. Four temperature sensors 
were mounted on the rotating apparatus, four others were attached to stationary components 
of the apparatus, and two sensors monitored the lab air temperature. Constant-temperature 
water from one Neslab RTE-221 unit flowed through copper pipes soldered to copper heat 
shields that surround the apparatus, and water from a second unit flowed through radiators 
with fans to hold constant the temperature of the room that houses the apparatus. The 
resulting thermal environment for the pendulum is constant to within a few mK per day. 

We consider two thermal effects. The first is the direct effect of a temperature change 
on the Iuj signal of the pendulum. We measured that a 1 K change in temperature of the 
apparatus caused a 9 ± 4 nrad change in the lu signal ((3 = 1.7 x 10~ 20 eV/K averaged 
over the sensors). Second, the temperature sensors on the rotating apparatus show small 
variations (~ 0.1 mK) at the turntable rotation frequency. Part of this signal may be due 
to friction in the turntable bearings which, if changing in time, could lead to changes in the 
lu signal of the pendulum. 

1. Astronomical signals 

We determined the direct temperature effect by fitting the measured mean temperatures 
of the sensors on the apparatus to the astronomical basis functions. We then multiplied 
these by the measured feedthrough, 1.7 x 10~ 20 eV/K, to obtain the systematic errors given 
as "temp, drift" in Table II. Similarly, we found the lu temperature effect by fitting 
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the lou signals of the rotating temperature sensors to the astronomical basis functions and 
multiplied by the feedthrough. We determined the feedthrough by performing a correlation 
analysis between the leu signals on the rotating temperature sensors and the leu signal of 
the pendulum. The feedthrough was found to be < 160 x 10 -20 eV/K. The systematic error 
is given as u lcu temp." in Table II. 

2. Lab-fixed signals 

The standard deviation of the average temperature of the sensors on the apparatus for 
the lab-fixed data runs was 2.8 mK. Because there was no correlation between the average 
temperature and the orientation of the spin pendulum within the apparatus, the direct tem- 
perature feedthrough created a random variation of the Iuj pendulum signal. The systematic 
error for the temperature drift effect, given in Table IV is 2.8 mKx 1.7 x 10~ 20 eV /K / V N — 1 
where N is the number of data runs. In measurements 1 and 2, N — 16, while measurement 3 
had N = 8. The standard deviation of leu signals on the rotating temperature sensors for the 
lab fixed data runs was 28 fiK. Again, these signals were not correlated with the orientation 
of the spin pendulum and would produce a random error. We used the same feedthrough 
for the bearing friction temperature effect as above for the astronomical signals to assign a 
systematic error for the la; temperature feedthrough of 28 /iK x 160 x 10~ 20 eV/K/^N — 1 
for the lab-fixed measurements. 

F. Spin content of the pendulum 

The uncertainty in the pendulum's spin content adds a scale-factor uncertainty to all our 
constraints. 

1. Lab- fixed signals 

Our constraints on A z , Czx, Czy, Czz and the CP- violating interaction of Eq. 8 with 
A << 1 A.U. are derived from laboratory-fixed signals. In this case we adopt the spin content 
and its 18% uncertainty given in Eq. 13. 
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2. Astronomical signals 



Our constraints on these signals are based on effects occuring at sidereal or solar fre- 
quencies, and steady lab-fixed /3's were ignored in the analysis. This allowed us to eval- 
uate the constraints using a more precise experimental value for the spin content based 
on the fact that it is extremely unlikely that new fundamental physics could have pro- 
duced the lab-fixed signal in Table V. This would require either a preferred-frame A 
that happened to point exactly opposite to the Earth's spin axis (this has a probability 
An/(4vr) w [5(3 E cos^/(2(3 N )] 2 = 7 x 1(T 5 ) or a CP-violating interaction of Eq. 8 with 
-Re < ^ < 3 x 10 6 km. If we neglect these unlikely possibilities, our measured torque 
/t(5 N = —(0.2537 ± 0.0057) fNm provides, via Eq. 34, an absolute calibration of the spin 
content of our pendulum, 



which gives a 2.3% scale-factor uncertainty for the astronomically signals. This was always 
negligible incomparison to the statistical error. 

VII. RESULTS 

In this section all constraints given in numerical form are la and include both random 
and systematic errors. Constraints presented in graphical form are at the 95% confidence 
level. 

A. Preferred-frame constraints 

1. Phenomenological cosmic preferred frames 

To constrain the parameters A and C we require results from lab-fixed as well as astro- 
nomically modulated signals because the Az, Czx, Czy and Czz signatures have no sidereal 
modulation. The lab fixed results are displayed in Table V; the combined result of the three 
measurements is 



AT^p = (9.80 ± 0.22) x 10 : 



,22 



(40) 



fa = (-0.03 ±0.28) x 10~ 20 eV 
fa = (0.00 ± 0.04) x 10~ 20 eV , 



(41) 
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where /3n = Ptss — Pn™, and the overall scale factor error from the uncertainty in N p is now 
included. The results in Table V allow us to extract values for either A z , or else for C Z x, 
Czy and Czz- The result in Eq. 41 implies that 

A z = /3 N /cos^ 

= (-0.4 ± 4.4) x KT 21 eV . (42) 

On the other hand, to constrain the boost terms we must invert the equations 

v x (ti) C zx +v Y (ti) CzY+vzih) C zz = /9 N (ii)/cos* 
vx(t 2 ) C Z x+v Y (t 2 ) C Z Y+vz(t 2 ) Czz = /3 N (t 2 )/cos^ 
vx{h) C Z x+v Y {h) CzY+v z {h) C Z z = Pn^)/ cos^ . 

where v is the laboratory velocity in units of c and U refers to the mean time of a lab-fixed 
measurement given in Table V. We obtain 

C zx = (-4.9 ± 8.8 ± 0.9) x 10~ 18 eV 
C ZY = (-9.6 ± 18.9 ± 8.8) x 10" 18 eV 

C zz = (+24.8 ± 38.6 ± 20.9) x 10" 18 eV , (43) 

where the second error reflects the scale factor uncertainty. 

Our complete constraint on the electron's Lorentz-symmetry violating rotation parameter 
A is displayed in Table VI; it is roughly two orders of magnitude more restrictive than 
previous work by Hou et al. [22], and substantially below the benchmark value mg/M Planck = 
2 x 10~ 17 eV. 

Our ler constraint on the electron's helicity-generating parameter B is 

B = (+0.50 ± 1.13) x 10~ 19 eV . (44) 

Table VII displays our limits on the electron's generalized helicity tensor C. We are not 
aware of any comparable measurements of either B or C. 

2. Standard Model Extension 

The minimal Standard Model Extension or SME incorporates the possibility of Lorentz 
and CPT violation by invoking pervasive, feeble, static spin 1 and spin 2 fields that are 
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TABLE V: Three separate measurements of lab-fixed signals. Units of torque, k5 and (3 = n5/N p 
are 10 -16 Nm and 10 -20 eV, respectively. Dates refer to the mean time of the measurements. 
Errors in (3 include the systematic and random errors given in Table IV but not the scale-factor 
uncertainty. 



signal 


15/10/06 value 




15/6/07 value 


8/3/08 value 


/?n 

(3 N - (3^° 


-2.49 ±0.11 
-1.62 ±0.07 
-1.62 

±0.00 ± 0.07 




-2.57 ±0.08 
-1.67 ±0.05 
-1.62 

-0.05 ±0.05 


-2.52 ±0.12 
-1.64 ±0.08 
-1.62 

-0.02 ±0.08 


(3e 


±0.04 ±0.11 
±0.03 ± 0.07 




-0.06 ±0.08 
-0.04 ±0.05 


±0.08 ±0.12 
±0.05 ±0.08 


TABLE VI: la constraints from our work and from Hou et al.[22] on the Lorentz- violating rotation 
parameters defined in Eq. 1. Units are 10~ 22 eV. The errors in Az are larger than those in Ax 
and Ay because of the greater systematic uncertainty in lab-fixed signals. It is assumed that C 
terms can be neglected. 


parameter 




this work 




Hou et al.[22] 


A x 
A Y 
A z 




-0.67 ± 1.31 
-0.18 ±1.32 
-4 ±44 




-108 ± 112 
-5± 156 
107 ± 2610 



defined in the heliocentric frame (for electrons in Minkowski spacetime 44 free parameters 
are involved). 

We follow conventional experimental practice [22-25] and first ignore possible violations 
of boost symmetry. Then our constraints on A (Table VI) translate directly to constraints 
on the SME parameters b e . Table VIII compares the electron's b parameters to the corre- 
sponding proton and neutron parameters determined from hydrogen maser[23, 24] and dual 
129 Xe/ 3 He maser measurements [25]. The neutron and electron parameters have comparable 
upper limits. The hydrogen maser results have an uncertainty 10 4 times greater so, in effect, 
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TABLE VII: la constraints on the C parameters of Eq. 3. The last three rows in this table 
are inferred from lab-fixed measurements. Units are 10~ 18 eV and the scale factor uncertainty is 
included. 

parameter value parameter value 



C X x +0.96 ±2.16 

C X y +1.84 ±4.24 

C X z -3.92 + 10.16 

C ZX -4.92 ± 8.86 

Czy -9.6 + 20.9 

C Z z +24.8 + 43.9 



C Y x -3.74 + 2.17 

C YY -7.76 ± 4.25 

C YZ +16.17 + 10.23 



they place a substantially weaker limit on b of the proton. 

TABLE VIII: la constraints on the Lorentz-violating rotation parameters for electrons, protons 
and neutrons. The SME boost terms are assumed to be negligible. Units are 10~ 22 eV. Proton 
and neutron results are taken from Refs. [23], [24] and Ref. [25], respectively. 

parameter electron proton neutron 



b x -0.67+1.31 < 2 x 10 4 0.22 + 0.79 

by -0.18 + 1.32 < 2 x 10 4 0.80 + 0.95 

b z -4 ± 44 



We constrained boost terms in the SME by reevaluating our limits after setting the 
velocity of the sun to zero. The resulting boost-dependent coefficients are denoted by primes. 
The sidereally modulated data were fitted with 6 free parameters, A X) A Y) C' xx , C' XY , C YX 
and C'yy- I n addition, we used the lab- fixed results in Table V to constrain A Z) C' zx and 

C' ZY = C' ZY + UneC' zz . (45) 

The A and C coefficients determine [26] the following combinations of SME parameters 

A x = b e x A Y = b e Y 
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C'xx = [br + d- ~ 2g c - 3g T + 4rf+ - d Q ]/2 
C'xy = {dxY — Hzt) ± tan c Hyt 
C'yx = Hzt 

C' yy = [2g c -g T -b T + 4rf+ - cL - d Q ]/2 + 

tane [d YZ - H X t\ 
C'zx = Hyt — dzx C' ZY = —H X t 

C' zz = [g T - 2d + + d Q ] . (46) 
Our values for these parameters are listed in Table IX. 

TABLE IX: la constraints on electron SME parameters from this experiment. Scale factor errors 
are included. The A coefficients listed here differ from those in Table VI because here we included 
additional boost terms whose functional forms are not fully orthogonal to those of the A coefficients. 

parameter electron value (GeV) 



A x 


(-0.9 ± 1.4) 


x 10- 31 


A Y 


(+0.9 ± 1.4) 


x 10" 31 


A z 


(-0.3 ±4.4) 


x 10 -30 




(+0.9 ±2.2) 


x 10- 27 


C'xy 


(+0.1 ± 1.8) 


x 10- 27 


C'yx 


(-4.1 ±2.4) 


x 10- 27 


C'yy 


(-0.8 ±2.0) 


x 10~ 27 


C'zx 


(-4.9 ±8.9) 


x 10- 27 


C'zY 


(+1.1 + 9.2) 


x 10- 27 



3. Ghost condensate gravity 

The result in Eq. 44 provides a limit on —M 2 /F (see Eq. 4) of 

M 2 



(-0.50 ± 1.10) x lCT 19 eV , (47) 
F 
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which improves on the value given in Ref. [6] by three orders of magnitude. The corre- 
sponding 95% confidence constraint on the "ghost condensate parameter" M/F is shown as 
a function of M in Fig. 13. A condensate-mediated spin-spin interaction with a strength 
comparable to gravity would have the ghost condensate parameter M/F ~ 10~ 16 [6]. 




M [evl 

FIG. 13: (color online) Upper limits on the "ghost condensate" parameter M/F as a function of 
M. The shaded area is excluded at 95% confidence. The vertical dashed lines show values of M 
having particular cosmological significance. 



4- Non- commutative geometry 

Inserting the constraint on \A\ from Table VI into Eq. 7, we obtain an la upper limit on 
|0|, the minimum observable patch of area, 

\&xz\ or \&yz\ < 4.9 x 1(T 59 m 2 [1 TeV/A] 2 (48) 
\&xy\ < 1.5 x 1CT 57 m 2 [1 TeV/A] 2 (49) 

which, assuming A = 1 TeV, corresponds to energy scales of 2.8 x 10 13 GeV for Qxz or Qyz 
and 5.0 x 10 12 GeV for Q XY . 

5. Axial torsion 

Many theorists (see Refs. [30] for a recent review) have speculated that the Reimann 
geometry of general relativity should be replaced by the more general Reimann-Cartan 
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geometry that contains an additional non-symmetric field called torsion. Torsion does not 
directly couple to unpolarized matter and hence is difficult to detect; however torsion's 
fully-antisymmetric axial component A is minimally coupled through covariant derivatives 
to fermion spins. Lammerzahl[31] noted that results from Hughes-Drever experiments can 
place useful bounds on A, and Kostelecky, Russell and Tasson[32] recently showed how 
experiments, such as the one we report here, determine all three spatial components of A. 
Following Ref. [32], we find for minimal coupling (Ref. [32] considers more complex cases 
as well) that our results in Table VI constrain any background torsion field (assumed to be 
roughly constant over the region of the solar system) to 

A x = -4A x /(3hc) = (+4.5 ± 8.9) x 1(T 16 nT 1 (50) 
A Y = -AA Y /(3hc) = (+1.2 ±8.9) x 1(T 16 m" 1 
A z = -4:A z /(3hc) = (+2.7 + 29.8) x KT 15 m" 1 . 

B. Boson-exchange constraints 

1. CP-violating monopole-dipole interactions 

The constraints on CP-violating monopole-dipole interactions extracted from our lab- 
fixed and solar and lunar-source data are shown in Fig. 14 (see also Table X). This new 
work improves over previous work[27-29] by factors of up to 10 4 . 

2. Velocity- dependent interactions 

The limits on velocity-dependent boson-exchange forces, extracted from our experimental 
bounds on interactions between the spin pendulum and the sun are given in Table X. Figure 
15 shows range-dependence of the velocity-dependent couplings f v and f±. We are not aware 
of any previous measurements of such forces. 

C. A test of the Equivalence Principle for spin? 

It may be of some interest to interpret our constraints as an "equivalence-principle" test 
for intrinsic spin, i.e. to ask if the gravitational mass of an electron with its spin pointing 
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TABLE X: la boson-exchange constraints from interactions with the Sun and Moon. Note that 
f v = 2</^<7y. The solar and lunar constraints assume A >> 1.5 x 10 11 m and A >> 4 x 10 8 m, 
respectively. 

parameter solar constraint lunar constraint 



&gg/(hc) (-3.5 ± 8.5) x 

f x / (he) (-0.1 ± 2.1) x 

9a9v /(fc) (+0.2 + 1.2) x 



10" 37 (+0.2 ± 1.6) x 10~ 34 

10~ 32 (-1.1 ±8.6) x 10" 29 

10" 56 (-3.1 ±2.4) x 10" 50 




10" 2 10"' 10° 10' 10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 11 10 12 
I [ml 



FIG. 14: (color online) Upper limits on \g^,g^\/{%c) as a function of interaction range A; the 
shaded region is excluded at 95% confidence. Our results and previous work by Youdin et al.[27], 
Ni et al.[28] and Wineland et al.[29] are indicated by solid, dash-dotted, dotted and dashed lines, 
respectively. Our work does not provide constraints for 10 km < A < 10 3 km because integration 
over the terrestrial surrounding is not reliable in this regime (see Ref. [20]). 

toward the sun is identical to that of an electron whose spin points away from the sun. In 
this case the relevant parameter is 

- _ l m g]t ~ [ m g]a _ 2/3 sun 

v ~ l(Wt + K]a) v N ' 1 } 

where the subscripts t and a refer to an electron whose spin points toward or away from the 
sun, (3 snn is the component of our signal (3 that tracks the sun, and Vn = —5.04 x 10~ 3 eV 
is the Newtonian potential energy of a laboratory electron in the field of the sun. We find 
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FIG. 15: (color online) Constraints on the velocity-dependent couplings f v and f± inferred from 
the infinite-range results in Table X. The shaded area is excluded at 95% confidence. In making 
this plot from the A = oo results in Table X we neglected the eccentricities of the earth's and 
moon's orbits. 

that /3 sun = (-0.62 ± 1.30) x 10" 22 eV, which leads to 

fj = (+2.4 ±5.2) x 10~ 20 , (52) 

where the errors are la. A considerably tighter limit comes from considering electrons falling 
toward the center of our galaxy. In this case Vjv = ^ e ©o where Go = 220 ± 20 km/s is the 
velocity of the solar system around the center of the galaxy, which gives Vjv = — 0.275±0.050 
eV. In this case we find (3 gal = (+0.27 ± 1.50) x 10" 22 eV, which leads to 

fj = (+0.2 + 1.1) x 10~ 21 . (53) 
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D. Measurement of the spin density in S111C05 



The absolute "gyrocompass" calibration of the spin content of our pendulum in Eq. 40 
provides a clean comparison of the spin density, p s , of room-temperature grade 22 SmCo 5 
to that of Alnico 5, 

p s (SmCo 5 ) - p s (Alnico) = -N^/(V V ) 

= (-3.66 ± 0.08) x 10 22 spins/cm 3 , (54) 

when the SmCo 5 and Alnico have identical magnetizations of B = 9.6 ± 0.2 kG. The 
quantities rj and V are defined in Eq. 10. We have neglected the very small measured 
magnetic moment (/Ib = —0.04) of Sm in SmCo 5 compared to the Co moment (— 8.97/xg), so 
that the equality of the magnetic fields in the two materials implies that the magnetizations 
satisfied M(Co) = M(Alnico) . The spin density in Alnico is 

p s (Alnico) = B /(fi fi B ))f Alnico 

= (7.85 ± 0.17) x 10 22 spins/cm 3 , (55) 

where /Alnico is taken from Eq. A2. Therefore 

p s (SmCo 5 ) = (4.19 ± 0.19) x 10 22 spins/cm 3 . (56) 

VIII. SUMMARY 

We have shown that a torsion balance fitted with a spin pendulum can achieve a constraint 
of ~ 10~ 22 eV on the energy required to flip an electron spin about an arbitrary direction fixed 
in inertial space. This is comparable to the electrostatic energy of two electrons separated 
by 90 AU. We then use this and related constraints to set sensitive limits on preferred- 
frame effects and non-commutative spacetime geometries. Our upper limits on the electron's 
rotation-noninvariant parameters \A e x Y \ < 1-5 x 10~ 22 eV and \A e z \ < 44 x 10~ 22 eV are 
substantially smaller than previous work[22] and up to 5 orders of magnitude below the 
benchmark value m 2 /M P i anck = 2 x 10~ 17 eV. Corresponding constraints on preferred-frame 
effects involving protons and neutrons are given in Refs. [23] and [25] . Interpreting our results 
as a constraint on non-commutative geometries we find that the minimum "observable" area 
is \®xz\ or \Qy,z\ < 4.9 x 10~ 59 m 2 , which corresponds to a length scale £ = 355 / G ut where 
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^gut = hc/(10 16 GeV). (Our limit \Qxy\ < 1-5 x 10 57 m 2 is weaker because it is derived 
from a lab- fixed signal which has a larger scale factor uncertainty.) These limits assume 
that the electron remains point-like up to an energy scale of 1 TeV. 

By analyzing 17 months of data, we obtained the first constraint on all 9 components of 
the Lorentz violating "boost" parameter C for any particle. Our upper limits on the absolute 
values of these components ranges between 2 x 1CT 18 and 23 x 1CT 18 eV. To facilitate SME 
analyses, we also present constraints on C", the tensor helicity with respect to the sun. 
Our electron constraints have roughly the same sensitivity as related neutron constraints 
obtained with a dual-gas maser[25]. 

Our constraints on CP-violating monopole-dipole interactions, derived from laboratory- 
fixed and solar-source signals, improve on previous work by factors of up to 10 4 , while our 
results for interactions between the spin pendulum and the sun or moon provide the first 
sensitive test for exotic spin-and-velocity-dependent interactions of electrons. 

Finally, our observation of the gyrocompass effect on the electron spins provides a precise 
result bearing on the electronic spin structure of Sm C05. If we assume that effects from 
exotic physics are negligible, our measurements determine the density of polarized electrons 
in Sm Co 5 is (4.19 ± 0.19) x 1(T 22 cm" 3 at a field of 9.6 kG. 

We are now taking a new series of measurements that probe the interaction between 
two spins. This work is motivated by the recent development of a consistent model of 
gravity with spontaneous Lorentz-symmetry violation [5] that predicts dramatic new spin- 
spin forces. These new experiments also probe the velocity-independent spin-spin forces [11] 
that can arise from generic boson exchange. 
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TABLE XI: Magneto-mechanical factors of Fe-Co alloys from Ref. [34]. Percentages of Fe and Co 
are by weight. 



%Fe %Co g' f 



75 25 1.918 ±0.002 0.957 

50 50 1.916 ±0.002 0.956 

25 75 1.902 ±0.002 0.949 



Astrophysics. CEC is grateful for an NSF Fellowship. 

APPENDIX A: SPIN CONTRIBUTIONS TO THE MAGNETIZATION IN 
SmCo 5 AND Alnico 5 

1. Alnico 5 

The composition of Alnico 5 by weight is 51% Fe, 24% Co, 14% Ni, 8% Al and 3% Cu and 
its magnetization arises primarily from the polarized spins of 3d electrons in the Fe and Co. 
The orbital moment of these electrons is small because of quenching in the inhomogeneous 
crystalline electric fields[33]. The polarized spins are confined to small, needle-shaped regions 
of Fe-Co alloy that precipitate out during heat treatment. Although we are not aware of 
direct measurements of the spin and orbital moments of Alnico 5, we can infer the spin 
fraction from "Einstein-De Haas" measurements [34] of the magneto-mechanical factors, g', 
of magnetically "softer" Fe-Co alloys. The g' measurements determine, /, the fractional 
spin contribution to the magnetic moment via 

/ = ^- (AD 

Table XI gives the measured spin fractions in various Fe-Co alloys. We adopt the mean and 
standard deviation of these values for the Alnico spin fraction, 

/Alnico = 0.953 ± 0.005 (A2) 
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TABLE XII: Magnetic moments and spin fractions for RC05 compounds. Unless otherwise noted 
results are deduced from neutron scattering measurements[35, 36]. 

compound A*(Coi) /(Coi) M c °n) /(Con) 



YCo 5 a 
CeCo 5 6 
SmCo 5 c 
Nd Co 5 c 



-1.77(2) n B 

-1.86(2) fi B 
-1.95(3) fi B 



0.74(5) 

0.76 
0.72(6) d 
0.67(5) 



-1.72(2) fiB 

-1.75(2) hb 
-1.90(3) hb 



0.84(4) 

0.90 
0.85(6) d 
0.80(3) 



a Ref. [35, 36] 

^NMR measurement from Rcf. [37] 
c Rcf. [36] 

d value inferred from neighboring compounds 



2. SmCo 5 

a. Co in Sm C05 

Although the magnetization of cobalt in SmCo5, as that of Alnico, arises primarily from 
polarized 3d electrons, the orbital moments are not as fully quenched because the rare earth 
atom alters the crystal lattice structure. Neutron[15, 35, 36] and photon[17] scattering 
studies, as well as NMR experiments [3 7], have shown that the Co atoms in RC05 (R is one 
of the rare earth elements) occupy two different crystal sites; 2 atoms in the 2c site (Coi) 
and 3 in the 3g site (Co n ). Although the spin fraction of the Co magnetization in SmCo 5 
has not been measured, experiments consistently find that the Co magnetic moments in 
RC05 compounds are independent of temperature and, as shown in Table XII, vary only 
slightly across the rare earth series. We assume that the Co spin fractions also vary only 
slightly with R to infer /(Coi) and /(Con) for Sm Co 5 as given in Table XII. The total spin 
moment of the 5 Co atoms 

^s(Co) = 2/(Co I )//(Co I )+3/(Co„)/x(Co„) 

= (-7.14 ± 0.39) f/ B , 
//(Co) = 2/i(Coi) + 3/i(Co n ) 
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which implies that 



= (-8.97 + 0.10) // B , 
fco = Ais(Co)MCo) = 0.80 ± 0.04 



(A3) 
(A4) 



b. Sm in Sm Co§ 

The Sm atoms in SmCo 5 are in the 3 + ionic state with a (4/) 5 (6s) 2 electronic config- 
uration. Hund's rules predict that the free Sm 3+ ion's orbital momentum has the largest 
possible antisymmetric value of L = 5 and that the spins are in the symmetric S = 5/2 state. 
Therefore one expects large orbital as well as spin contributions to the magnetic moment. 
One also expects that the ion's ground state has S and L coupled to the minimum J = 5/2. 
This simple model predicts // = —0.714 // B , fx s = +3.57 /ib and //£ = —4.29 //b- However, 
the Sm ions are embedded in a lattice and the crystalline and exchange fields alter this 
simple picture. Neutron scattering experiments [15] on SmCo 5 find that //(Sm) = — 0.38 //b 
at a temperature of 4.2K. However, at room temperature, ksT becomes comparable to the 
splitting of the magnetic substates and thermal population of the various Mj states causes 
the measured Sm magnetic moment to be very small, //(Sm) = — 0.04/i#[15]. Using the 
Sm wavefunctions in Ref. [15], we calculate that at the temperature maintained during our 
experiments 



The X-ray magnetic Compton scattering measurement reported in Ref. [17] probes the 
spin moment of Sm in Sm C05 directly, but provides no information about the total magnetic 
moment. The ratio of Co to Sm spin moments from this measurement, R = —0.23 ± 0.04, 
together with Eq. A3, corresponds to 



c. Ratio of Sm to Co spin moments 

Our calculation of the number of polarized electrons in our pendulum, Eq. 13, depends 
on R, the ratio of Sm to Co spin moments. The neutron scattering results in Eqs. A5 and 
A3 give 



/is(Sm) = +3.56//B 



(A5) 



/i S (Sm) = +1.67 ±0.31 fi B . 



(A6) 



R = // s (Sm)/// s (Co) = -0.50 . 



(A7) 
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On the other hand the photon scattering work[17] directly gives 

R = -0.23 ± 0.04 . (A8) 

We accomodate this difference by assuming that R is equally likely to lie anywhere between 
these two values. 

APPENDIX B: "FEETBACK" TILT ELIMINATION SYSTEM 

We developed a digital feedback system that continuously corrects for the varying tilt of 
the laboratory floor and imperfections in the turntable itself. It analyses inputs from an 
orthogonal pair of co-rotating electronic tilt sensors and, every 2 complete revolutions of the 
turntable, computes the tilt as a function of turntable angle and cancels this predicted tilt 
by controlling the length of feet that support the turntable. The performance of the leveling 
system is shown in Fig. 16 which demonstrates that the system compensates not only for 
the slowly-varying tilt of the laboratory floor but also for the more rapid variations caused 
by imperfections in the turntable bearing. 

The tilt is measured by Applied Geomechanics Inc. inclinonometers (AGI's) mounted on 
the rotating apparatus close to the upper attachment point of the pendulum's suspension 
fiber. A constant tilt of the rotation axis appears as a lu signal on each AGI, with a 90° 
phase difference between the two units. The analog signals of the AGI's are digitized by the 
data acquisition system and analyzed to determine the tilt with a precision of ~ 1 /irad/ v^Hz. 
This tilt is then removed by adjusting the lengths of 2 of the 3 feet upon which the apparatus 
rests. The lengths of the two feet are are controlled by varying their temperatures. Figure 17 
shows a cross section through one foot. The expanding and shrinking components of a foot 
consist of two lead rings (1) which are soldered to a copper disk (2). Thermal energy can be 
pumped into or out of the copper disk by a Peltier element (3) thermally coupled to a brass 
block (4) that is held at constant temperature by circulating water from a temperature- 
stabilized reservoir. Two G10 rings (5) thermally isolate the lead rings from the laboratory 
floor and from a stainless steel disk on top, upon which the apparatus rests. The Peltier 
element and the brass block are clamped to the copper disk by one bolt (6). Specially formed 
G10 pieces provide thermal isolation between the bolt, the copper disk and the brass block. 
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FIG. 16: Performance of the "feetback system" showing the tilt measured by co-rotating level 
sensors. Upper panel: "feetback" off. The rapid fluctuations arise from imperfections in the 
turntable bearing, slower variations arise from the varying tilt of the laboratory floor. Lower 
panel: "feetback" on. 

The thermal capacity and thermal resistance of the copper disk and the lead rings produce 
a 14-second delay in the expansion of the lead rings compared to the change in the heat flux 
provided by the Peltier element. In addition, the thermal conductivity of the Peltier element, 
which connects the copper disk to the temperature reservoir, affects the time constant of 
the response. Taking these factors into account, we developed a model, shown in Fig. 18, 
to predict a foot's response to a change in heat flux delivered by its Peltier element. The 
model poses no problem for a real-time calculation because it can be solved analytically. 
The delayed response of the feet together with an 8-second integration time of the AGI's 
leads to a low-pass behavior of the leveling system. Therefore the response of the feet to fast 
changes of the tilt caused by imperfections of the ball bearings in the turntable is limited. 
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FIG. 17: Cross section of a foot of the Eot-Wash rotating torsion balance. 

The tilt caused by the imperfections of the bearing is periodic in the turntable angle. This 
allowed us to express the bearing tilt in Fourier coefficients of the turntable angle. Once 
these Fourier coefficients are calculated, they are fed forward to compensate for the bearing 
imperfections. A simplified flow chart of the feedback loop is shown in Fig. 18. With the 
feetback switched on, a typical tilt of a few /zrad per day is reduced to less than 10 nrad. 

After the simple model described above was implemented, we found it necessary to add 
one further refinement. The DC offsets of the AGI's drift slowly in time. A feedback system 
that holds constant the AGI signals does not guarantee true level. With the feedback on, a 
drifting offset produces a lu AGI signal with the opposite phase lag between the two sensors 
compared to a true tilt. Our model extracts this Fourier component and uses it to determine 
the AGI readings that correspond to true level. 



52 



simple thermal model 




heat flux peltier element 
heat capacity copper disk 
heat capacity lead rings 
heat resistance peltier element 
heat reststance copper disk 



turntable angle 



raw 
AGI^ 

raw 
AGI^ 



calibration 



AGI 



calibration 



transformation 
into lab frame 



projection 
on foot A 
and foot B 



At 



T1A . 



T2A , 



hfA 



ATA 




AGI Y 



AGIX 



AGI Y 



calculation 
of fourier 
coefficients 
correction 



solve simple 
thermal model 
for foot A 



,T1A 



,T2A 



2 <■ 

ir. I ^» 



hfA 



calculate feed 
forward signal 
from fourier 
coefficients 




update fourier 
coefficient at 
first zero crossing 
after 2 turns 



to solve model for foot B . 



inverse 

peltier element 
function 



to peltier 
element 



DAC calibraton 
|update DACl — >• supply 



FIG. 18: Thermal model and flow chart of the simplified feetback loop implemented at the Eot- 
Wash II rotating torsion balance 
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